Resonant States in the Electronic Structure of the High Performance Thermoelectrics 
AgPb^SbTe2+m ; The Role of Ag-Sb Microstructures 
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Ah initio electronic structure calculations based on gradient corrected density functional theory 
were performed on a class of novel quaternary compounds AgPbmSbTe2+m, which were found to 
be excellent high temperature thermoelctrics with large figure of merit ZT^2.2 at 800K. We find 
that resonant states appear near the top of the valence and bottom of the conduction bands of bulk 
PbTe when Ag and Sb replace Pb. These states can be understood in terms of modified Te-Ag(Sb) 
bonds. Electronic structure near the gap depends sensitively on the microstructural arrangements 
of Ag-Sb atoms, suggesting that large ZT values may originate from the nature of these ordering 
arrangements. 
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Thermoelectrics (TE) drive wide variety of applica- 
tions related to solid state refrigeration and small scale 
power generation. Over last four decades narrow band- 
gap semiconductor alloys based on Bi-Te compounds for 
cooling, and Sii-a^Ge^c and Pb-Te alloys for power gen- 
eration have dominated technological applications. In 
recent years there has been a strong impetus to improve 
the efficiency of TE and this has led to a constant search 
of new materials and new concepts. 

Success in discovering novel TE hinges on the ability 
to achieve the rather challenging task of synthesizing ma- 
terials with simultaneously high electronic conductivity 
(cr), high thermopower (5) and low thermal conductiv- 
ity {hz). These solid state properties define the figure of 
merit ZT=(j^ T/a^, where T is the operating tempera- 
ture. Several new ideas have been proposed to achieve 
high values (>1). One is the so-called electron crys- 
tal phonon glass (ECPG) where one reduces the lattice 
thermal conductivity without dramatically reducing the 
power factor (a^) Q. The other idea is to increase 
crS^ by manipulating electronic density of states (DOS) 
using 0, 1, and 2 dimensional quantum confinement ef- 
fects I3. In general aS^ depends on the DOS through 
the transport distribution function Si(e) given by: 



(1) 



where i is the transport direction, the summation is over 
the first Brillouin zone, Vi{k) is the velocity of the carriers 
with wave vector /c, and r{k) is the relaxation time. Yet 
another concept to increase ZT by introducing a narrow 
T^i{e) centered near the chemical potential was suggested 
by Mahan and Sofo 3] . The last two concepts introduce 
resonant-like states near the Fermi energy Ej. 

Recently a new family of complex chalcogenide com- 
pounds, AgPb^SbTe2+m was described; several mem- 
bers of this family (m=18) when doped appropriately 
exhibit large values (-2.2) at 800K jj]. These com- 



pounds were originally designed to have an average cubic 
NaCl structure {Fm3m symmetry) with Ag, Pb and Sb 
atoms being statistically disordered on the Na sites. In 
fact powder X-ray diffraction studies support this model. 
However our careful single crystal X-ray diffraction and 
electron diffraction studies reveal that is not the case (see 
Fig.^). There is clear experimental evidence that the Ag 
and Sb atoms are not statistically disordered with the Pb 
atoms. Instead there is a strong driving force that causes 
long range ordering in the crystal, the nature of which 
depends on the value of m and the exact experimental 
conditions with which the materials were prepared l5|. 
In fact it is possible to create a statistical disorder of all 
cations on the Pb sites in the samples (such as quenching 
from a melt) but such samples have inferior TE proper- 
ties. 




FIG. 1: Electron diffraction pattern showing a lowering of the 
crystal symmetry from cubic to orthorhombic. It corresponds 
to [112] zone axis in fee symmetry with cubic unit cell (a=6.4 
A). Superstructure spots can be indexed in orthorhombic 
symmetry according to ao ^9.1 A, bo ^15.8 A and co ^^22. 3 
A unit cell. The corresponding zone axis is [012] 0. Indexation 
is given for both symmetries. 

Given the lower crystallographic symmetry of 
AgPbmSbTe2+m we cannot regard these systems as solid 
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solutions between AgSbTe2 and PbTe but as bona-fide 
quaternary compounds. This is not surprising since 
there is a strong enthalpic driving force to not disperse 
Ag+ and Sb^+ randomly in the structure lest we create 
coulombic instabilities as we have argued earlier 0j. In 
addition, it is not correct to view the Ag and Sb atoms as 
dopants because they are present in large stoichiometric 
proportions. The AgPbmSbTe2+m materials achieve 
a higher ZT than PbTe at elevated temperatures by 
exhibiting different temperature dependence in all three 
properties influencing ZT 4]. Therefore there is a 
substantial modulating effect on the TE properties of 
PbTe when the Ag and Sb atoms adopt the proper 
ordering patterns in it. In this letter we discuss the elec- 
tronic structures of several AgPbmSbTe2+m compounds 
(m=10, 16, 18, and 30) based on different ordering 
models to investigate this issue. The results of our calcu- 
lations demonstrate that significant qualitative changes 
occur in the electronic DOS near Ej in AgPbmSbTe2+m 
vis-a-vis PbTe that could cause favorable enhancements 
in <jS^ and its temperature dependence. 

Given the close structural relationship of 
AgPb^SbTe2+m to PbTe, a fundamental question 
is how the PbTe electronic structure gets modified by 
extensive substitution with Ag and Sb atoms and their 
microstructural arrangements. In order to address this 
we have performed first ah initio electronic structure cal- 
culations using full-potential density-functional method 
for a series of different microstructural arrangements. 
We have used the WIEN2K package for our calculations 
1^. This is a linearized full-potential augmented plus 
local orbitals method within the density functional 
theory formalism 0. We took the generalized gra- 
dient approximation in Ref. |^ for the exchange and 
correlation potential. Scalar relativistic corrections were 
included and spin-orbit interaction was incorporated 
using a second variational procedure |^ . Convergence of 
the self-consistent iterations was performed using 10 (for 
single and Ag-Sb pair atoms'''), 15 (for layer of Ag-Sb 
pairs), 18 (for chain of Ag-Sb pairs) and 2 (for cluster of 
Ag-Sb pairs) k points inside the reduced Brillouin zones 
to within 0.0001 Ry with a cutoff of -6.0 Ry between the 
valence and the core states. 

First, we performed calculations with isolated Ag and 
Sb atoms in the PbTe lattice in order to obtain a clear 
picture of their individual role in modifying the elec- 
tronic structure of PbTe (Fig. [2ti). Then both Ag and 
Sb were introduced to simulate stoichiometrics relevant 
to those of AgPbmSbTe2+m compounds. Given that the 
exact crystal structure is not known, several plausible 
microstructural models were examined, all of which in- 
volved long range ordering of the atoms. In one model the 
Ag and Sb atoms were placed in monolayers (Fig. [33, c). 
In another model the Ag and Sb atoms were placed along 
straight infinite chains running parallel to a crystallo- 
graphic unit cell axis (e.g. c-axis)(Fig. I21i). In yet a 
third arrangement, the atoms were placed in the center 
of a 3x3x3 supercell to create a "AgSbTe2 nanodot" em- 
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FIG. 2: Unit cell models for (a) single Ag atom 
AgPb3iTe32 , (b) Ag-Sb layer perpendicular to [001] direc- 
tion in AgSbPbi8Te2o, (c) Ag-Sb layer perpendicular to the 
fee [111] direction in AgSbPbioTei2, (d) Ag-Sb chain parallel 
to the [001] direction in AgSbPb3oTe32, and (e) Ag-Sb cluster 
in AgSbPbieTeis. For the reason of clarity we show only Pb 
fee lattices with Pb in blue, Ag in red, and Sb in green colors. 



bedded in a PbTe matrix (Fig.I2t)- The chain and cluster 
models are in qualitative agreement with our experimen- 
tal transmission electron microscopy (TEM) observations 
of Ag-Sb ordering 4]. Although we believe these ar- 
rangements capture much of the crystal physics in these 
materials, we recognize that there are many more that 
could be considered which nevertheless should result in 
similar general conclusions. 

To model the Ag(Sb) isolated and Ag-Sb pair atoms 
we constructed 2x2x2 supercells with 64 atoms. For the 
isolated case, we chose Ag(Sb) at the origin of the super- 
cell with a separation of two lattice constants (12.924 A) 
between the Ag(Sb) atoms (Fig.[2t)- For the Ag-Sb pair, 
we considered two arrangements (not shown in Fig. [21), 
one where the Ag and Sb are far apart (Sb at the origin 
and Ag at the center of the supercell) with a separation 
distance of ~11.19 A and the other where the Ag and Sb 
are as close as possible (Sb at the origin and Ag at the 
next nearest neighbour site of Sb) with a separation dis- 
tance of ~4.57 A. For the structure where Ag-Sb layers 
are separated by several Pb layers we also considered two 
cases where the Ag-Sb layer is normal to the [001] direc- 
tion in a 1x1x5 supercell (40 atoms/cell) with the Ag-Sb 
layer located in the z=0 plane (Fig. [^d) and where the 
Ag-Sb layer is normal to the [111] direction (Fig. [2t). 
The fee unit cell can be viewed along the [111] direction 
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as a hexagonal unit cell. We have used a 2x2x2 hexagonal 
supercell (48 atoms/cell) with the Ag-Sb layer perpendic- 
ular to the c-axis which is the [111] direction in the fee 
unit cell. To model the chains we used a 2x2x2 supercell 
where the Ag-Sb chains are oriented parallel to the [001] 
direction (Fig. Eli) and separated by 12.924 A. Finally, 
for the "AgSbTe2" clusters we constructed a 3x3x3 su- 
percell (216 atoms/cell) and the cluster consists of six 
Ag-Sb pairs located at the center of the supercell with 
a minimum separation distance between two clusters of 
-12.924 A (Fig.Et). 

The nature of defect states in narrow band gap semi- 
conductors in general and PbTe in particular are impor- 
tant problems but still not completely understood. Lent 
et. al. 10] argued that in PbTe, due to its large dielec- 
tric constant Coulomb effects associated with charged 
impurities were screened out and local bonding effects 
determined the nature of the defect states. Instead of 
donor or acceptor states appearing in the band gap as 
seen in conventional wide band gap semiconductors they 
argued that one should see resonance states outside of the 
band gap region. In this letter we provide for the first 
time theoretical justification of this idea starting from 
first principles. As discussed below, our ah initio calcu- 
lations clearly show that the states associated with iso- 
lated Ag(Sb) atoms in PbTe are not only resonant with 
the valence band (VB) and conduction band (CB) states 
of PbTe but they also reduce the band gap. 

The total DOS for isolated Ag atoms is shown in 
Fig. int- It can be seen that Ag introduces states near 
the top of the PbTe VB. Partial DOS analysis shows that 
these states consist mostly of p orbitals of the six nearest 
neighbor Te atoms of Ag. These states are resonant with 
the VB and extend into the PbTe gap region. On the 
other hand the isolated Sb single atoms introduce reso- 
nant states near the bottom of the PbTe CB (Fig. ^p) 
which extend nearly —0.75 eV into the CB starting from 
the CB bottom. Sb and its Te nearest neighbours atoms 
have the highest contribution to these resonant states. 
The Sb p states hybridize with Te p states in the range 
(-0.25,0.5) eV. Therefore these states are not only reso- 
nant with the PbTe CB but they extend into the PbTe 
gap. 

Results for isolated Ag-Sb pairs are consistent with 
the Ag(Sb) single atom results in the sense that Ag in- 
troduces new states near the top of VB, whereas Sb in- 
troduces new states near the bottom of CB (Fig. |nfc,d) 
decreasing the PbTe gap. Both cases (Ag-Sb far apart 
and Ag-Sb next nearest neighbours) (not shown in Fig.[2|) 
show semiconducting behaviour with a very small gap 
and a more rapidly increasing DOS near the VB and CB 
extrema as compared to the DOS of PbTe. The spe- 
cific features of the DOS in the gap region are very dif- 
ferent for these two cases. Total energy comparison for 
AgSbPb3oTe32 shows that these two structures are very 
close in energy, the case when Ag-Sb atoms are far apart 
has a lower energy by —20 me V/ (unit cell). 

It is interesting to compare the DOS results for the 
layer structures of Ag-Sb (Fig. [2l3,c). When the Ag-Sb 



layer is perpendicular to the [001] direction, the states 
associated with the Ag-Sb layer completely fill the PbTe 
gap giving a semimetallic behaviour (Fig. ^), whereas 
the Ag-Sb layer perpendicular to [111] direction show a 
semiconducting behaviour (Fig.^t*)- This indicates that 
the electronic structure of AgPbmSbTe2+m systems and 
consequently the electronic properties are very sensitive 
to the microstructural arrangements of Ag-Sb atoms. 
The Ag-Sb chain model shows semiconducting behaviour 
(Fig. ^). This chain model has the same stoichiometry 
(AgSbPb3oTe32) as the Ag-Sb pair models. Total energy 
comparisons show that the chain model has a lower 
energy by 0.2 eV/(unit cell) than the Ag-Sb pair models 
suggesting that Ag-Sb chain orderings along [001] direc- 
tions are favorable microstructures. This is consistent 
with the results of electron crystallographic studies 
which indicate the presence of Ag-Sb chains in the 
crystal 1^]. The DOS results for the "AgSbTe2" cluster 
model show also semiconductor behaviour (Fig.^). 
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FIG. 3: Total DOS of (a) single Ag atom, (b) single Sb atom, 
(c) Ag-Sb pair with the Ag-Sb distance of —11.19 A, and (d) 
Ag-Sb pair with Ag-Sb distance of —4.57 A, in PbTe. For 
comparison the total DOS of PbTe is shown in dashed line. 
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FIG. 4: Total DOS of (a) Ag-Sb layer model perpendicular to [001] direction, (b) Ag-Sb layer model perpendicular to [111] 
direction, (c) Ag-Sb chain model, and (d) Ag-Sb cluster model. The total DOS of PbTe is shown in dashed line. 



The calculations for different microstructural arrange- 
ments of Ag-Sb in PbTe show a generic feature; when 
Sb atoms replace Pb atoms, Sb hybridize with Te atoms 
forming strong covalent interactions. When Ag replace 
Pb, the p states of Te which are the nearest neighbours 
of Ag are strongly perturbed. Therefore the electronic 
structure of AgPbrnSbTe2+m compounds depend sensi- 
tively of these perturbed Te states by the Ag(Sb) atoms. 

We have shown for the first time that the details of 
the DOS near the energy gap of AgPbmSbTe2+m depend 
sensitively on the microstructural ordering arrangements 
of Ag-Sb pairs in PbTe. The nature of the states near the 
top VB and bottom CB in these quaternary compounds 



is substantially different than in PbTe. The common 
feature of these Ag-Sb arrangements is that they have a 
more rapidly increasing DOS near the gap as compared 
to bulk PbTe due to the appearance of distinct resonant 
states. It is well accepted that resonant structures in the 
DOS near E/, created by quantum size effects B [ill. 
superlattice engineering |l3 or chemical means [T^ 
are very desirable features because they could enhance 
the TE figure of merit ZT. 
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